Diffusion-limited exciton-exciton annihilation in single-walled carbon nanotubes: 

A time-dependent analysis 
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To provide physical insight into the recently observed photoluminescence saturation behaviors in 
single-walled carbon nanotubes implying the existence of an upper limit of exciton densities, we 
have performed a time-dependent theoretical study of diffusion-limited exciton-exciton annihilation 
in the general context of reaction-diffusion processes, for which exact treatments exist. By including 
the radiative recombination decay as a Poissonian process in the exactly-solvable problem of one- 
dimensional diffusion-driven two-particle annihilation, we were able to correctly model the dynamics 
of excitons as a function of time with different initial densities, which in turn allowed us to reproduce 
the experimentally observed photoluminescence saturation behavior at high exciton densities. We 
also performed Monte Carlo simulations of the purely stochastic, Brownian diffusive motion of one- 
dimensional excitons, which validated our analytical results. Finally, we consider the temperature- 
dependence of this diffusion-limited exciton-exciton annihilation and point out that high excitonic 
densities in SWNTs could be achieved at low temperature in an external magnetic field. 

PACS numbers: 78.67.Ch,71.35.-y,78.55.-m 



I. INTRODUCTION 

Excitons in single-walled carbon nanotubes (SWNTs) 
are stable quasi-particles with large binding energies and 
significantly influence their interband optical propertiesii 
However, they have been reported to be rather efficiently 
eliminated at high densities through the exciton-exciton 
annihilation (EE A) process;^ although their emission and 
absorption energies remain stable even at high densities."^ 
Recently, the intensity of photoluminescence (PL) from 
SWNTs was found to saturate at high pump fluence,— 
implying the existence of an upper limit in the density 
of excitons, which was estimated to be an order of mag- 
nitude smaller than the expected Mott density. The ex- 
istence of such an upper limit, which poses a significant 
hinderance in the observation of lasing, a Mott transition, 
or excitonic Bose-Einstein condensation in SWNTs, was 
attributed to efficient EEA facilitated by the diffusive 
motion of excitonsi^i^ii 

The dynamics of diffusion-limited EEA can be ana- 
lyzed in the general context of reaction-diffusion pro- 
cesses, which have been extensively studied by physi- 
cists, chemists, biologists, and ecologists and serve as 
simple models for studying a variety of non-equilibrium 
problems . Moreover, it has been shown that 
such simple diffusion-driven reactions exhibit interest- 
ing non-equilibrium phase transitions and universality 
classes,"i2. with connections to many-body theory.— In 
such models, particles, or "agents," of one or more species 
execute random walk in d-dimensions (where d could also 
be fractional in the case of fractal geometry) and undergo 
reactions upon collisions, leading to changes in their pop- 
ulation often accompanied by an appearance and disap- 
pearance of various phases. Such systems often exhibit 
rich phase diagrams that can be fully studied with nu- 
merical simulations even when analytical solutions are 
not available. For example, a widely studied reaction is 



the two-particle annihilation given hy A + A ^ 0, where 
the interaction is assumed to be of the "hard-core" type, 
leading to the mutual destruction of two particles upon 
collision.^ Starting with an initial population Ng, one can 
analyze the ensemble averaged population Nt at a given 
time t. Clearly, this is a non-equilibrium many-particle 
process that is driven by noise, and its only steady state 
is achieved when the population vanishes. 

An interesting feature of diffusion-driven reactions is 
the presence of spatio-temporal fluctuations in such pro- 
cesses that cannot be ignored especially at lower dimen- 
sions and lead to the breakdown of mean-fleld type as- 
sumptions. For example, diffusion-driven two-particle 
annihilation can be written in differential form as follows: 



dtinixA)) = DV^n{x,t)) - {n^x,t)), 



(1) 



where (.) stands for an ensemble average. The first term 
on the right denotes the diffusion process, while the sec- 
ond term is for two-particle annihilation and thus has 
quadratic dependence. Note that the annihilation term 
has an average of n'^{x, t) whereas mean- field theory will 
simplify this to (n(a;,t))^, thus neglecting fluctuations 
of the form {n'^{x,i)) — (n(a;,t))^. In one dimension. 



the population asymptotically decays as t 



-1/2 



power- 



law whereas the mean-fleld theory, which ignores fluctu- 
ations, predicts a faster decay of This discrepancy 
originates from the fact that d-dimensional diffusion with 
d < 2 is recurrent and the particles return to their pre- 
vious position with high probability^lS Hence, the reac- 
tion is slowed down leading to a smaller exponent in the 
power-law decay. Indeed, the mean-field result is recov- 
ered in three dimensions, which is above the critical di- 
mension dc — 2 for this problemiii Thus, an exact treat- 
ment of even such a simple process requires the inclusion 
of correlations. 

In this paper, we undertake a time-dependent study 
of the diffusion-limited EEA process in the presence of 
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radiative decay. In particular, we consider the following 
two coupled and competing reaction-diffusion processes 
in f-D to model the dynamics of excitons in nanotubes 
at various densities: 

A + A^kA (fc = 0,l); (2a) 
A^B (2b) 

where A represents excitons and B photons. The first 
equation represents exciton-exciton annihilation, which 
is either complete [k — 0) or partial [k — I), while the 
second reaction is just the radiative decay of excitons 
with radiative lifetime = f/7r- It is noteworthy that 
only the first reaction is diffusion-driven whereas as the 
radiative decay takes place independently. This diffu- 
sive motion of excitons is due to the random collisions 
with phonons^ In this sense, the annihilation reaction is 
driven by diffusive noise, which we assume to be of the 
Gaussian form, whereas the radiative decay is governed 
by a Poissonian noise and, hence, is a pure jump process. 
We consider a simple diffusion process in which the diffu- 
sion constant D is independent of the spatial and tempo- 
ral coordinates. Furthermore, the excitonic dimension is 
assumed to be much smaller than the nanotube dimen- 
sion. We are interested in determining the population 
of both species as a function of time. The population of 
species S, or photons, is proportional to the PL intensity 
measured experimentally. In particular, we wish to know 
the fraction of population which decays radiatively and 
how this fraction changes as the initial population is in- 
creased. As one can imagine, upon increasing the initial 
density of excitons in the 1-D nanotube the annihilation 
reaction becomes more efficient whereas the radiative de- 
cay rate can be safely assumed to be independent of the 
density. As we show later, this leads to saturation of the 
PL intensity as the initial population density is increased, 
which is consistent with the experimental observations.^ 



II. SOLUTION USING FIRST-PASSAGE 
DISTRIBUTION 

The case of two-particle annihilation without the ra- 
diative decay has been extensively studied, and exact 
results for the population as a function of time are 
known Fi2ii2ii2£i2i Here we use the first-passage time dis- 
tribution of Brownian motion to first study the anni- 
hilation reaction without decay. We recover the exact 
analytical results for this case before proceeding to in- 
clude the radiative decay term. We begin by deriving 
the exact result for a single pair, or the "independent 
pairs" case, and use it to obtain an approximate solution 
for the many-particle, or the "correlated" case. Monte 
Carlo simulations are performed to check the validity of 
our results. This purely stochastic method employing 
the first-passage time distribution is a simple and natu- 
ral way to study the annihilation reaction as the collisions 
which drive the reactions must obey such distributions. 



A. Annihilation without decay; k = 

Consider A'o pairs of species A randomly arranged 
on a line of length L and executing Brownian diffusion 
with diffusion constant D. Let us first consider just the 
two-particle annihilation process without the decay as in 
Eq. ((2^) with k = Q. Let nyi(t) denote the average frac- 
tion of initial population which is still "alive" at time t. 
We keep the discussion in this section as general as pos- 
sible without explicitly identifying species A ov B unless 
absolutely required. 

In 1-D, only the nearest neighbors at any given time 
can undergo annihilation due to restrictions placed by 
lower dimensionality. This prompts us to first consider 
the case for a single pair of particles and generalize the 
result to the many-particle case. As this is equivalent to 
different pairs annihilating independently of each other, 
we refer to it as the "independent pair" case. Let c?o be 
the initial distance between the pair and P(c?o, t) denote 
the survival probability of this pair at time t. To cal- 
culate P{do,t), we need to find the probability that a 
pair with initial distance do does not undergo collision 
till time t. As both particles are executing independent 
Brownian motion, their relative motion is also Brown- 
ian with diffusion constant 2D which starts at do. Thus, 
we need to find the probability that a Brownian motion 
starting at do does not reach zero till time t. This can be 
readily found from the first-passage time distribution of 
a Brownian motion as^^ 

where erf(.) is the error function. For the independent 
pair case, the fraction of population that is alive at time 
t, nA{t), reads 

n^(i) = Vcrf(^)=Verf(^), (4) 

where the distance between the particles of the i-th pair is 
do,i at time t = 0. We have imposed a periodic boundary 
condition making the line into a ring without any loss of 
generality. As A^o tends to infinity, the above sum can 
be expressed as an integral over the distribution of do.iS. 
We restrict ourselves to the case when the particles are 
randomly arranged on the line at i = 0. Thus, do,iS which 
are the nearest neighbor distances can be thought of as 
the "waiting times" for a Poisson process and have an 
exponential distribution with mean do = L/2Nq. 

For the many-particle case we realize that there are 
twice as many ways for a pair to annihilate as for the in- 
dependent case due to the presence of two nearest neigh- 
bors for each particle, and hence, the mean distance for 
the correlated case is just a half of the independent case 
in the large Nq limit. For this limit, the exact result can 
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be obtained as 



due to mutual independence. Thus, 



dxdexp (—3x) erf 
^ ^ \2VDtJ 



X 



exp {P^Dt) erfc ( VfPDt 



(5) 



which is the result of Torney et al}^ with f3 — 4No/L. In 
the asymptotic limit, Eq. (O yields a power-law decay of 
as mentioned earlier. 
At long times the initial correlations between the par- 
ticles are completely wiped out, and this power-law de- 
cay is expected, irrespective of the initial distribution of 
particles. Such a power-law behavior has been recently 
observed in time-resolved transient absorption measure- 
ments on SWNTs^ 



B. Annihilation with decay; k = 0: Independent 
pair case 

Next, we include radiative decay of Eq. ([2)3) in our 
model and compute the population fraction of species A 
and i? as a function of time. As before, we first derive 
the exact result for the case of a single pair and use it 
as a kernel to express the result for independent pairs 
uniformly distributed along the tube. For a single pair 
separated by a distance do at t — 0, the surviving popu- 
lation fraction at t can be simply written as 



(6) 



where t) denotes the probability of n surviving par- 
ticles at time t, which remains to be calculated. Let us 
compute p{2,t), which is the probability that both par- 
ticles comprising the pair are alive at time t. Such a 
case is possible only if neither particle undergoes radia- 
tive decay or collision till time t. As the radiative decay 
of particles occurs independently of one another and also 
of the diffusion driven collision, we can simply multiply 
the individual probabilities to get 



p(2, t) = exp {—2-frt) erf 



do 



2VDt 



(7) 



Recall that radiative decay is a Poisson process with pa- 
rameter 7j. and the probability of it not happening till 
time t is exp(— 7^^). To compute p{l,t), which is the 
probability that exactly one particle out of the pair sur- 
vives, we realize that such a scenario is possible only if 
there is exactly one radiative decay in the time interval 
[0,t], say at t = r, and no collision before r. As in the 
interval (r, t] collisions cannot take place due to an in- 
sufficient number of particles for the reaction, we only 
include the probability of collision not taking place be- 
fore r. The probability that either of the particles decay 
in an infinitesimal interval dr about t is 2"frdT. As be- 
fore, we can multiply the probabilities for each sub-event 



do 



- / exp (— 27rT) erf 
Jo \2VDt 

X {2jrdT) exp {—2jr {t — t)) 



(8) 



From Eqs. ©, Q, and ([H]), nA{t) for a single pair can 
be obtained. As before, for the case of iVo independent 
pairs, we average nA{t) over an exponential distribution 
with mean do = L/2No. After some straightforward but 
tedious algebra we obtain 



nA{t) 



exp (-7,.i) 



exp 



1 



7,.t j erfc [y-frt/v 



-l-Vi/erf (VTrt) — v\ , 



(9) 



where we have introduced a dimensionless parameter v 
= Tu/Tr with td being the "diffusional time" d^/D. We 
emphasize that this is an exact result for the case of in- 
dependently colliding pairs that also undergo radiative 
decay. The above equation is valid only when v < 1, or, 
in other words, when radiative decay is slower. In the 
limit of extremely dilute initial population density, no 
annihilation can take place, and only radiative decay oc- 
curs. nA{t) would be simply exp(— 7ri) in that case. In 
the opposite limit when radiative decay rate vanishes, 
Eq. ^ indeed recovers the result of Eq. ([5]) , as expected. 

Figure [T] compares the result of Eq. ([9]) with Monte 
Carlo simulations, done by simulating the Brownian dif- 
fusive motion of each independent pair of particles, vali- 
dating our results. As u is the only physical parameter in 
the problem, changing do and D but keeping u constant 
should not alter the result, which was indeed confirmed 
by simulations. This fact should remain true even for the 
case of correlated pairs. 

Let us consider the behavior of nA{t) in the long and 
short time limits. For i ^ 1, only the radiative decay 
should dominate as the density of particle becomes too 
low to participate in annihilation. Thus, an exponential 
decay is expected. Taking limits explicitly, one obtains 



nA{t 



exp (-7r) 
1 + 



(10) 



As V is less than unity, so is the intercept of the above ex- 
ponential decay, hinting at the superexponential decay at 
short times due to annihilation. In the short time limit, 
when V < \ one expects only annihilation to dominate 
the decay, and one gets 



nA{t ^ 0) = 1 - 2yJ"frt/TTV. 



(11) 



which is indeed faster than an exponential decay as t — > 
0. Indeed, Fig. [1] confirms these findings. As the time 
progresses, the density of particles decreases due to de- 
creasing population, which slows down the annihilation 
reaction as it strongly depends on the density of the par- 
ticles. The radiative decay rate, on the other hand, is 
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FIG. 1: (color online). Comparison of Monte Carlo simula- 
tions (solid line) and exact analytical result (dashed line) for 
the fc = 0, "independent" case with radiative decay [Eq. (|9])]. 
The fraction of A [in (a)] and B [in (b)] populations is plotted 
as a function of dimensionless time Dt/do. The value of D = 
100, Tr — 80 and do — 20 = 0.05) were used for illustrative 
purposes. 



fixed, and thus, a crossover from annihilation dominated 
decay to a purely exponential radiative decay is expected. 
It can be defined to take place when dt = L/2NonA{t) 
becomes equal to v. This time t* is implicitly given as 



(12) 



As the initial density is increased, r* decreases finally 
vanishes at very high density, implying a purely expo- 
nential decay at all times. 

In order to calculate nB{t), we note that at any given 
time the rate of radiative decay is proportional to the 
instantaneous population of A, nA{t). In other words. 



(13) 



Hence, nsit) can be obtained from Eq. ([9]), upon direct 
integration, as 



nsit) = 7r / dm Air). 
Jo 



(14) 



In particular, the fraction of total population that decays 
radiatively is given by 



nsioo) = 



1 



1 + l/y^ 



(15) 



In Fig. [21 we use Eq. to plot the total number den- 
sity of species or photons, created as a function of 
initial density of species A, or excitons. A saturation be- 
havior of the number of photons created as the density 
of excitons is increased is seen even for the independent 
pair model agreeing with our intuitive understanding and 
experimental observations. Thus, the independent pair 
model captures all the essential features of the process. 
More importantly, the study of independent pair model 



0.1 0.2 0.3 

Initial Density A (a.u.) 



0.4 



FIG. 2: (color online). Population per unit length of B plot- 
ted against the initial density of A using Eq. (|15p to show the 
saturation behavior in PL intensity predicted by our model. 



identifies the relevant parameters of the process and the 
scaling relationships they must obey. It also provides the 
decay regimes that are relevant for each type of reaction 
viz., annihilation and radiative decay. The use of purely 
stochastic first-passage time distribution makes the so- 
lution transparent and simple, relying on the properties 
of diffusion rather than other formal methods, which al- 
though more general are less intuitive. 



C. Annihilation with decay; fc = 0: 

case 



Correlated pair 



As for the case of no decay, we scale the mean separa- 
tion between the particles by a factor of two in order to 
obtain a solution for the correlated case. This approxi- 
mate solution and the Monte Carlo simulations for the 
correlated case are compared in Fig. [3] The approximate 
solution agrees well with the simulations. A possible rea- 
son for the slower decay of analytical result compared to 
the exact result could be the following: For the single 
pair case, if one of the particle decays before undergoing 
collision, the remaining particle must decay radiatively 
and cannot undergo annihilation. However, for the cor- 
related case, this is not true as long as there are other 
neighboring particles and so annihilation becomes possi- 
ble. 



D. Annihilation with/without decay; k = 1 

The case of partial annihilation [fc = 1 in Eq. ^ii)] can 
be understood in terms of the results for fc = 0. Both 
processes are completely identical besides the fact that 
the annihilation in fc = 1 is half as slow as the fc = 
case. Consequently, if the initial density for the partial 
annihilation case is twice as much as the complete anni- 
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FIG. 3: (color online). Comparison of Monte Carlo simula- 
tions (solid line) and approximate analytical result (dashed 
line) for fc = "correlated pair" case with radiative decay 
(Eq. |9]). (a) The population fraction of A (a) and B (b) as a 
function of "dimensionless" time. Approximate result agrees 
fairly well with the simulations. The parameters for simu- 
lations are the same as in Fig. [1] except with do = 10 = 
0.025). Some reasons for disagreement with the simulations 
are discussed in the text. 
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FIG. 4; (color online). Comparison of fc = 1 (solid line) and k 
= (dashed line) for different values of do for the "correlated" 
case with radiative decay, (a) The population fraction of A 
(a) and B (b) as a function of "dimensionless" time, do values 
for k = 1 case simulations are half of the k = case shown on 
the graph). Rest of the parameters are the same as in Fig. [1] 



hilation case, one expects the two decays to be identical. 
Thus, the result for the k = 1 case can be obtained from 
Eq. ([5]) by replacing do with do/2. Even in the presence 
of radiative decay, the above argument should be true as 
the radiative decay occurs completely independently of 
the annihilation reaction. We verify this heuristic rea- 
soning by Monte Carlo simulations, as shown in Fig. 3] 



III. EFFECT OF TEMPERATURE ON THE 
UPPER LIMIT OF EXCITON DENSITY 

Thus far, we have not included the effect of tempera- 
ture (T) in our discussions. The temperature-dependence 
of the diffusion constant D can be approximated through 
the Einstein relation. 



D 



Mr' 



(16) 



where M is the mass of the exciton and F is the exciton- 
phonon scattering rate. At low temperatures, acoustic 
phonon scattering is expected to be dominant, and F for 
1-D is then given in terms of the deformation potential 
Ddr, as 



F 



(17) 



where p is the mass density and Vg is the sound velocity 
in the material. Thus, the temperature dependence of D 
due to acoustic phonon scattering is D oc a/T/A/'^/^. As 
the temperature is increased, F would become a sum of 
both acoustic and optical phonon scattering rates. 

As the relevant quantity in our model is v and not 
just -D, we need to include the temperature dependence 
of Tr in order to fully understand the temperature de- 
pendence of the EEA process in the presence of radiative 
recombination. The radiative lifetime Tr of a single 1-D 
exicton band is predicted to scale as^^ \/T, leading to v 
oc 1/T or the exciton diffusion length Ix oc VT. Thus, 
the saturation of photoluminescence due to EEA would 
become less effective at lower temperatures. However, in 
the case of carbon nanotubes, the presence of optically in- 
active, or "dark," states lying below the optically active, 
or "bright," state causes radiative lifetime to increase 
at low temperatures,— and this could favor the EEA 
process depending on the exact temperature dependence. 
By applying symmetry breaking perturbations such as a 
magnetic field, the dark state can be brightened;^ restor- 
ing the temperature dependence for the case of a single 
1-D exciton band at higher fields. Under such conditions, 
it may be possible to attain the Mott density of exci- 
tons in carbon nanotubes. Sustaining such high densities 
of excitons is the first step for any lasing applications 
and for observing excitonic Bose-Einstein condensation 
in carbon nanotubes. In addition, at lower temperatures 
exciton localization due to impurity traps or defects could 
completely stop the diffusive motion of excitons, further 
enabling the attainment of the Mott density.?! 

Finally, in our model, we have assumed a completely 
random motion of excitons, which leads to ordinary dif- 
fusion based on a Gaussian kernel. This assumption can 
also break down at lower temperatures or in other sce- 
narios when D becomes position- or density-dependent, 
leading to anomalous diffusion and changing the time- 
dependence of excitonic population. A time-resolved ex- 
periment, probing the excitonic or the photon population 
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at different temperatures, exciton densities, and mag- 
netic fields, can not only verify the validity of this model 
but also provide further insight into the EEA process in 
carbon nanotubes. 



IV. SUMMARY 

To provide fundamental physical insights into the re- 
cently observed photoluminescence saturation behaviors 
in single- walled carbon nanotubes, we studied the dif- 
fusion and two-particle annihilation of one-dimensional 
excitons in the general context of reaction-diffusion pro- 
cesses, for which exact treatments exist. By including 
the radiative recombination decay as a Poissonian pro- 
cess in the exactly-solvable problem of one-dimensional 
diffusion- driven two-particle annihilation, we were able to 



correctly simulate the density of excitons in single- walled 
carbon nanotubes as a function of time and density. 
Monte Carlo simulations were also performed by simu- 
lating the purely stochastic, Brownian diffusive motion 
of one-dimensional excitons, validating our results. Fi- 
nally, we discussed the temperature dependence of EEA 
and proposed possible experiments to verify the validity 
of this model. 
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